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ABSTRACT 

We present a method for solving the hghtcurve of an echpsing binary system 
which contains a Ccphcid variable as one of its components as well as the solutions 
for three eclipsing Cepheids in the Large Magellanic Cloud (LMC). A geometric 
model is constructed in which the component stars are assumed to be spheri- 
cal and on circular orbits. The emergent system flux is computed as a function 
of time, with the intrinsic variations in temperature and radius of the Cepheid 
treated self-consistently. Fitting the adopted model to photometric observations, 
incorporating data from multiple bandpasses, yields a single parameter set best 
describing the system. This method is apphed to three echpsing Cepheid systems 
from the MACHO Project LMC database: MACHO ID's 6.6454.5, 78.6338.24 
and 81.8997.87. A best-fit value is obtained for each system's orbital period and 
inclination and for the relative radius, color and limb-darkening coefficients of 
each star. Pulsation periods and parameterizations of the intrinsic color varia- 
tions of the Cepheids are also obtained and the amplitude of the radial pulsation 
of each Cepheid is measured directly. The system 6.6454.5 is found to contain 
a 4.97-day Cepheid, which cannot be definitely classified as Type I or Type 
H, with an unexpectedly brighter companion. The system 78.6338.24 consists 
of a 17.7-day, W Vir Class Type II Cepheid with a smaller, dimmer companion. 
The system 81.8997.87 contains an intermediate-mass, 2.03-day overtone Cepheid 
with a dimmer, red giant secondary. 

Subject headings: Magellanic Clouds — Cepheids — stars: AGB and post-AGB 
— stars: oscillations — binaries: eclipsing 

1. INTRODUCTION 

Large scale microlensing surveys have provided unprecedented resources for variable 
star research. Their long time baseline and stable, acciuate photometry are ideal for the 
detection and analysis of such objects and the large number of systems observed increases the 



Serena, Chile 

Email: nsuntzeff@noao.edu 

^°Department of Physics, University of Oxford, Oxford 0X1 3RH, U.K. 
Email: w . sutherlandOphysics . ox . ac . uk 



-4- 



probability of finding astrophysically-interesting objects that either have escaped detection 
or do not exist in our own galaxy. 

Very few regularly-pulsating stars are known to belong to eclipsing systems. One of 
the best studied is AB Gas, an Algol-type binary system, which contains a 5 Scuti type 
primary (Rodriguez et al. 1998). Rodriguez & Breger (2001) hst 6 additional S Scutis which 
are members of eclipsing binaries. The lone Galactic candidate for an eclipsing binary 
containing a Gepheid variable was BM Gas (Thicsscn 1956) but further study revealed that 
the variable was unlikely to be a Gcphcid (Fcrnic & Evans 1997). The astrophysical benefits 
of a Gepheid variable in an eclipsing system could be considerable. If the system is double- 
lined, a determination of the Gepheid's luminosity and mass can be made that is not only 
more accurate than existing measurements but also truly independent of the intervening 
steps in the distance ladder. Such a system would provide the most direct measurement of 
the mass of a Gepheid and would offer an independent calibration of the period-luminosity 
relation. 

Here we present the results of lightcurve analyses of three eclipsing Gepheid systems 
in the MAGHO project Large Magellanic Gloud (LMG) database: MAGHO ID's 6.6454.5, 
78.6338.24 and 81.8997.87. In Sec. 2 we describe the sources of the photometric observations. 
Sec. 3 describes the model used to generate the lightcurve of an eclipsing Gepheid system 
and Sec. 4 describes the inverse problem of computing the parameters from an observed 
lightcurve. In Sec. 5 we present the results obtained for the three systems. Finally, Sec. 6 
summarizes the analysis, describes work in progress and suggests future avenues for research. 

2. OBSERVATIONS 

Observations from several different sources are used in the analysis presented here but 
the majority are from the MAGHO Project database. The MAGHO observations were made 
with the refurbished 1.27m Great Melbourne Telescope at Mount Stromlo, near Ganberra, 
AGT, Australia. It has been equipped with a prime focus reimagcr-corrcctor with an integral 
dichroic beamsplitter to give a 0.5 sq. deg field of view in two passbands simultaneously: 
a 450-590 nm MAGHO V filter and a 590-780 nm MAGHO R filter. These are each sam- 
pled with a 2x2 array of 2048x2048 Loral GGDs which are read out concurrently via two 
amphfiers per GGD in about 70 seconds. The resulting images cover 0.5 square degrees 
with 0.63 arcsec pixels (Alcock ct al. 1995). Data reduction is performed automatically by 
Sodophot, a derivative of DoPhot (Schcchter, Mateo, & Saha 1993). MAGHO photometry is 
then transformed into Gousins V and R bands for further interpretation (Alcock ct al. 1999). 
The eclipsing, pulsating nature of systems 6.6454.5 and 78.6338.24 was identified through 
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visual inspection of their lightcurves. 

The echpsing nature of 81.8997.87 was first reported by the Optical Gravitational Lens- 
ing Experiment (OGLE) project (Udalski et al. 1999) (OGLE ID LMC_SC16 119952) and 
the hghtcurves for this system are a combination of observations in the OGLE and MACHO 
databases. OGLE observations were taken on the 1.3 m Warsaw telescope at Las Campanas 

Observatory, Chile operated by the Carnegie Institution of Washington. Photometry is in 
the standard BVI bands with the majority of the observations in the / band (Udalski et 
al. 1999). V and / measurements for 6.6454.5 (LMC_SC21 40876) have also been used here. 
The OGLE observations of the April 2000 eclipse of this system were graciously provided by 
Andrzej Udalski in advance of publication. 

For the 1999 eclipse, observations of 6.6454.5 were taken by Nick Suntzeff from March 
16 to April 21 on the 0.9m telescope at the Cerro Tololo Inter-American Observatory in 
Chile. The telescope is a 0.9-meter Ritchey-Chretien Cassegrain telescope with a dedi- 
cated 2048x2048 CCD detector. Observations were in standard BVI and were reduced by 
DAOPHOT, ALLSTAR and ALLFRAME. 

A complete listing of the observations, from all sources, used in this paper is available 
at http:/ /www. macho. mcmaster.ca/Data/EclCep/ecLcep. htm. For guidance regarding form 
and content a partial hsting is given in Tables 1-8. 



To analyse the lightcurve of eclipsing Ccphcid systems, a geometric model based on 
classical equations (Binnendijk 1960, for instance) was developed. In this model the stars 
are assumed to be traveling on circular orbits. This isn't the most general case but may 
be appropriate here. It is found that two out of the three systems studied here have their 
primary and secondary minima separated by almost exactly one half period, consistent with 
a circular orbit. The stars themselves are treated as circular disks with radius (R) and 
central surface brightness (Jo). A linear limb-darkening law of the form 



is adopted where x\ is the limb-darkening coefficient and 7 is the angle between the surface 
normal at that point and the line of sight. Integrating over 7 yields the light received from 
a non-eclipsed star as being proportional to 



3. 



MODEL 



J = Jo{l - xx + xx cos(7)). 



(1) 




(2) 
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As absolute dimensions of the stars cannot be extracted from the hghtcurve alone, the above 
parameters all represent ratios of observed quantities: the radii are measured in units of the 
orbital separation of the two stars, the surface brightnesses are normalized such that the 
mean light received outside of eclipse (Li + L2) is 1.0. 

In many programs the values of Jo are modelled directly with a different value for each 
bandpass. As we have simultaneous MACHO V and R observations for all systems at all 
phases, we fit the colour, specifically the ratio Jy / Jr for each star, rather than the surface 
brightness. From this surface brightness ratio the usual V — R colour index can be calculated. 
If I band observations arc available, as they arc for two of the three systems, the fiux ratio 
Jv/Ji must also be defined, but is not used as a fitted parameter. Instead the V — / colour 
index is computed from the V — R colour using a linear relationship derived from the VRI 
observations of Cousins (1980) of a sample of southern stars. Prom the colour index, Jy/ Ji 
can be calculated. To convert the surface brightness ratios into values for the individual 
surface brightnesses we employ the linear relationship between V surface brightness (Jy) 
and the V — R colour index identified by Barnes, Evans, & Parsons (1976) (hereafter B- 
E). This relationship holds for stars of all spectral types and for all pulsational phases of 
the Ccpheid. Unfortunately, it is not certain that the B-E relation applies to metal-poor 
Type II Cepheids and this is a potential source of systematic error in our results. However, 
not all Type II Cepheids are metal poor: one of the few Type II Cepheids known to be in 
spectroscopic binaries, the 2.4-day Type II Cepheid AU Peg (Harris, Olszewski, & Wallerstein 
1984), has an [Pe/H] = -1-0.1. The B-E relation apphes to Johnson filter definitions so the 
transformations of Bell & Gustafsson (1980) were used to convert to forms apphcable to our 
Cousins bandpasses. With values of Jy the other surface brightnesses follow from the ratios. 

The intrinsic variabihty of the Cepheids is handled by taking it! and Jv/Jr as the 
minimum radius and the mean colour, respectively. Variability is then added to these values 
in some functional form. The colour variation is parameterized by a third-order Pourier 
series: 



This produces a curve which, despite its discontinuity at Ai? = 0, reproduces the broad 
features of the radial variation curve. 




(3) 



The radial variation is modeled by the expression: 




(4) 



During eclipses, the decrease in the amount of light received from the system is calcu- 
lated analytically, based on the radii of two overlapping discs and geometrical considerations 
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(Binnendijk 1960). The apparent separation of the centres of the two stars is given by 

5^ — cos^ i + sin^ i sin^ 6 (5) 

Whether the system is in echpse or not is determined by comparing this distance to the 
instantaneous radii of the stars. Within echpse, the area obscured by the echpsing star is 
given by 

Area = -Rl (2q;i - sin(2Q;i)) + ]-Rl {2a2 - sin(2Q;2)) (6) 
/ \ Rf — Ro + S'^ 

, , Rr, — R^ -\- 5 , , 

= 2rIs 

Multiplying the area by the surface brightness of the eclipsed star gives the decrease in flux. 

The effects of limb- darkening within the eclipses are included by dividing the stellar disk 
into 100 concentric rings in a manner similar to Nelson & Davis (1972). A surface brightness 
is then assigned to each ring based on equation 1. The total light received from each star 
can then be calculated by computing the contribution from each ring and integrating over all 
of the rings. Within eclipses the integration must be performed within suitable limits based 
on the amount of each stellar surface that is visible. In practice, it is faster to compute the 
normal fluxes from equation 2 and then calculate the amount of fiux lost based on the eclipse 
geometry. This is done by using equations 6-8 repeatedly for eclipsed stars of varying radii, 
which correspond to the rings into which the stellar disk is divided, and then summing over 
the total number of rings. So the fiux lost would be given by 

100 

Ji {Area{ri) - Area{ri_i j) (9) 

1=0 

where ro to rioo are the radii of individual rings ranging from to the radius of the eclipsed 
star (Rl), Area{ri) is the area of a disk of radius covered by the eclipsing star, given by 
equations 6-8, and Jj is the surface brightness of that ring given by equation 1. This flux lost 
is then subtracted from the flux given by equation 2 to give the eclipsed star's contribution to 
the sytem flux. Integration over the eclipsed region would be more accurate in principle, but 
poses the practical problem of being more difficult to implement. The limits of integration 
are particularly difficult to express for a star of varying radius as Ri (as well as Jo in equation 
1) will be a function of time. Summation is more straightforward to implement and the loss 
of accuracy in the final flux value is, with the number of annuli being used, on the order of 
one part in 1000, as determined by comparison of trial summations over a full, uneclipsed 
disk with the result from the analytic expression (equation 2). 
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4. FITTING PROCEDURE 

Prior to the fitting of the adopted model to the observed hghtcurves, a prehminary 
analysis of each lightcurve was performed. This was done to produce a set of initial param- 
eters as a starting point for the fitting process. The orbital period can easily be estimated 
by measuring the average length of time between eclipses. Then, a coarse removal of the 
pulsation contribution to the lightcurve is achieved by phasing the data to the orbital period 
and fitting hyperbolas to both the primary and secondary echpses. The out-of-eclipse flux is 
assumed to be 1.0. The result is an averaged lightcurve which gives clearly deflned estimates 
for the average minima of both eclipses and average phase of external contact. Values for 
the inclination, radii and relative luminosities were obtained using an algorithm from Riazi 
(1992): from the values for the system light, in two filters, at the eclipse minima, this algo- 
rithm gives the relative luminosities of the two stars as well as the ratio of radii. From these, 
and the approximate phase angle at external contact, the inclination and relative radii can 
be computed. These values are then improved by performing a fit to minimize relative 
to a simplified model. Once the boundaries of the eclipses have been defined the parameters 
controlling the Cepheid variability can be estimated by examining only the out-of-eclipse 
data. The pulsation period is found using a routine period from Press et al. (1992) on the 
out-of-eclipse data points. This routine uses an algorithm by Lomb (1976) to compute the 
power present at various frequencies in unevenly spaced data. Once the pulsation period is 
known, the data outside of eclipse can be phased to it and a Fourier series fit to the flux 
variation. The initial radial variation is simply set to be one quarter of the radius of the 
variable star and the offset is set to be one quarter of a period later than the time zeropoint 
for the lightcurve. The parameter set resulting from the above procedures will not accurately 
describe the lightcurve but it is usually sufficient to produce a lightcurve close enough to 
that observed to allow fltting to proceed and to be in the same valley in parameter space 
as the global minimum. 

The data are then subjected to a minimization by means of a standard Levenberg- 
Marquandt method from Press et al. (1992). This procedure alternates between two com- 
plementary methods for flnding a minimum x^: the steepest descent method far from the 
minimum, switching smoothly to the inverse Hessian method as the minimum is approached. 
This approach combines the advantages of both methods: rapid progress is made towards a 
minimum and, once located, accurate determination of the best-flt parameters is achieved. 
The drawback is the possibility of termination within a local minimum rather than the true 
minimum. This possibility can be reduced through reasonable selection of the initial values 
for the fltting process and examination of the flnal parameter set. An array of flags allows 
each parameter to be flt or held flxed. The data from both fllters (or all three fllters, when 
available) are flt simultaneously, allowing parameters such as the inclination and the radii. 
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which correspond to physical dimensions of the system (and so should not vary from filter 
to filter) to be determined from all of the data. 

The limb-darkening coefficients prove to be the parameters most difficult to fit and have 
to be treated separately. A fit is first performed with all hmb-darkening coefficients held 
fixed at their initial values, usually 0.5, to produce improved values for the other parameters. 
Another fit is then performed with the limb-darkening coefficients allowed to vary to produce 
optimum values for all parameters. The x will not vary far from 0.5, typically staying within 
0.4 - 0.7, the range expected for x of most main-sequence and giant stars. The relative 
movement of the x\ values can then be examined to see if it matches expectations based on 
the relative temperatures and surface gravities of the two stars. In practice the uncertainties 
on the limb-darkening coefficients prove to be too large for meaningful analysis. More "in 
eclipse" data would improve the constraints on these coefficients. 

To distinguish between primary- variable and secondary- variable configurations, we solved 
for the system parameters with either configuration in turn. The set with the lower was 
adopted. 

We also investigated the possibility of sources of flux in each system in addition to the 
two eclipsing components. This involved assuming a third light source, of various bright- 
nesss and refitting the system. In each of the three systems this produced no significant 
improvement in the fit. 

5. RESULTS AND DISCUSSION 

The best fit parameters for each of the systems are shown in Tables 9-11. The orbital 
and pulsational periods are both in days and the inclination is given in degrees. The radii (R) 
and amplitude of the radial variation (ARamp) are relative to the orbital separation of the 
two stars. ARshift, which measures the shift of the radial change relative to the temperature 
change, has units of days. Only the ratio of surface brightnesses Jy/ Jr is tabulated as it was 
the only one that was fit directly and the other surface brightnesses are computed from it as 
outlined above. The limb-darkening coefficients, defined in equation 1. Once the 

relative surface brightness and radius of each star has been determined they can be combined 
with the mean system magnitude to compute the magnitude of each star in all filters. Also 
computed for each of the three Cepheids is the value of Wr = R — 4.0(V^ — R), an index 
which corrects for most of the effects of reddening and effective temperature differences. 

Figures 1-3 show the primary eclipses for each system along with the best fit lightcurve. 
Figure 4 shows a period-H^R diagram (P-L diagram) for MACHO Cepheids and the locations 
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of the three Cepheids studied here. This diagram is essentially free of reddening and allows 
us to classify the three Cepheids under study based on their relation to other LMC Cepheids 
without an explicit correction for extinction. The uncertainty ranges in the magnitude values 
in Tables 9-11 and Figure 4 were estimated from the range of possible component magnitudes 
based on the uncertainties in the best fit parameters. These are statistical uncertainties and 
likely underestimate the true uncertainties. 

The colors and magnitudes of each star allow some general comments on the evolutionary 
state of each system if we assume each star follows a standard, single-star evolutionary 
history. This assumption is not unwarranted here given the large orbital periods of the 
systems. First, a crude correction for extinction must be made. To account for foreground 
reddening values of E(B — V) are adopted from the map of Galactic foreground color excess 
toward the LMC pubhshed by Schwering & Israel (1991). This yields values of E{B — V) for 
the three systems as follows: 0.08 for 6.6454.5, 0.08 for 78.6338.24 and 0.10 for 81.8997.87. 
These allow us to determine values for Ay and Af{ when combined with the standard value 
of the ratio of total to selective extinction, Ry = Ay / E{B — V) = "i.l (Cousins 1980), and 
Ar/Av = 0.77 for our Cousins R and V from the interstellar extinction relations of Cardelli, 
Clayton, & Mathis (1989). 

This procedure is clearly sufficient for the system 78.6338.24, producing Vo = 16.20±0.03 
which agrees with the period- relation of Alcock ct al. (1998) for type 11 stars with log 
P > 1.1. Applying the correction described above to the system 81.8997.87 fails to produce 
V and R magnitudes for the Ccphcid that are consistent with those expected of an overtone 
Cepheid of its period. Given the system's proximity to the 30 Doradus star-forming region it 
is not unreasonable to expect substantial extinction along this line of sight within the LMC 
and the closest Cepheid on the same plate, the overtone 81.8997.128, also appears well below 
the overtone band in the P-L diagram prior to applying a reddening correction. The period- 
magnitude relation for overtone Cepheids of Baraffe & Alibert (2001) gives K = 15.86 
for the Ccphcid in 81.8997.87, implying Ay = 1.31, a total value of E{B - V) = 0.42 
and Ro = 15.56. Both Vo and Ro are consistent with the overtone bands in the period-V 
magnitude and period-R magnitude plots for MACHO LMC Cepheids. A total value of 
E{B — V)= 0.41 was obtained by de Marchi et al. (1993) for selected regions of 30 Doradus. 

A further correction may also be necessary for the system 6.6454.5, however, the am- 
biguity in the classification of the Ccphcid (see below) precludes definitively comparing its 
properties to those expected from a period-luminosity relation as was done with the other 
two systems. In light of this we adopt only the correction for foreground reddening. 

After correction for reddening the values of V and V — R are converted to the L — Teff 
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plane by assuming ijllmc = 18.5 mag and using 

logiTeff) = 4.199 - ^0-08369 + 0.3493(y - R) (10) 

a transformation of the Chiosi, Wood, & Capitanio (1993) semi-empirical calibration. Figures 
5 and 6 show theoretical isochrones representative of the metallicities of Type I and Type II 
stellar populations: Y=0.25, Z=0.008 isochrones from Bertelli et al. (1994) in Figure 5 and 
Y=0.230, Z=0.0004 isochrones from Fagotto, Bressan, Bertelli, & Chiosi (1994) in Figure 
6, along with the properties of the three systems. The error bars are computed exclusively 
from the errors in the magnitudes and colors given in Tables 9-11. 

Based on the properties tabulated in Table 9 the system 6.6454.5 is found to contain 
a Cepheid as the secondary with a brighter, bluer primary. The nature of the Cepheid is 
unclear as its location in the P-L diagram (Figure 4) places it between the Type 1 and 
Type 11 bands, inconsistent with both classifications. The large amplitude of the radial 
variation, 0.323 ± 0.008 of the Cepheid's minimum radius is more consistent with a Type II 
classification. Fundamental-mode Type I Cepheids have typical radial variations of of 10% 
or less (Armstrong et al. 2001) while Type II Cepheids show larger radial excursions in the 
range of 30-50% of the minimum radius (Lebre & Gillet 1992). 

If the Cepheid is assumed to be a Type II Cepheid it is either making an excursion 
from the AGB or moving off the HB to the AGB (a less hkely scenario given its period). 
The companion, displayed in Figure 6, which also appears to be considerably evolved, is too 
luminous to fit either of these scenarios. It is possible that this system is not in the LMC but 
is instead a foreground object. A reduction in the assumed distance to the system to 17.8 
kpc is necessary to shift the Cepheid's properties to fit a post-HB evolutionary state. By 
contrast, if the system is compared to Type 1 isochrones (Figure 5) the Cepheid's location is 
consistent with that expected but the companion appears to be too bhic to fit the isochrones. 
If the Cepheid is indeed Type 1, an additional reddening correction would need to be applied 
to make its location in Figure 4 consistent with the fundamental. Type I band. This would 
imply an even higher effective temperature and luminosity for the companion. 

The status of 78.6338.24 is less ambiguous. It consists of a Type II Cepheid secondary 
with a hotter, but somewhat dimmer primary. With a pulsational period in excess of 17.5 
days the Cepheid would classified as a W Vir type, which is consistent with its large radius 
relative to its companion and with its location in the P-L diagram (Figure 4). This system 
presented several challenges to modeUing. The pulsation period of the Cepheid was found to 
not be constant over the duration of the observations. Despite being very small in magnitude 
(less than 1% of the period) this drift in period produced a substantial decrease in the quality 
of the fit. It was corrected to some degree by assuming a pulsational angular frequency that 
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was a slowly varying fuction of time, parameterized by: 

CO = ujo + AiBt + A2iBtf (11) 

where the are parameters to be determined by fitting and B is a constant, set by trial 
and error, to ensure that the Ai are of the same order as the other parameters. Inspection 
of the complete set of residuals after fitting revealed indications of non-sphericity in one (or 
both) of the system components. Both the asphcricity and the "period drift" are consistent 
with the large radius and tenuous outer envelope of a W Vir star. Their locations in Figure 
6 show both components to be wcll-cvolved, post-HB or post-AGB objects. 

The binary 81.8997.87 is distinct from the other two systems in several ways. Its vari- 
able is an intermediate-mass Cepheid pulsating in the first-overtone mode. Furthermore, the 
Cepheid is the primary with a considerably cooler, dimmer companion. The amplitude of 
radial pulsations is low, 0.060 ± 0.006 of the minimum Cepheid radius, as expected for an 
overtone Cepheid. Gieren (1982) provides an exphcit determination of the radial displace- 
ment for the galactic overtone Cepheid SU Cas. The radial amphtude of that star is 0.026 
of the mean Cepheid radius, similar to the value we find for 81.8997.87 (note that workers 
tend to avoid Baade-Wesslink analyses of overtone Cepheids due to the small dynamic range 
of the obscrvables). Table 5a of Pel (1978) lists physical properties derived from fitting 
model atmospheres to the continuum colors of a sample of Galactic Cepheids. The values 
of AR/ (R) obtained for the six overtones in the sample range from 0.042 to 0.080. Once 
photometric contamination from the companion is removed the Cepheid's lightcurve shape, 
parameterized by the Fourier ratio i?2i, is consistent with those of other MACHO overtone 
Cepheids of similar period. Its location in the P-L diagram (Figure 4) further reinforces this 
view. The companion's properties suggest a RGB star, possibly K class or later, however 
this combination of binary components is in poor agreement with current models of single 
star evolution. This is reflected in the disagreement with the isochrones seen in Figure 5. 
Indeed any system consisting only of these two stars may be evolutionarily inconsistent. 

The observational coverage of this system is far from ideal. The 800-day orbital period 
limits the number of primary eclipses in the MACHO database to only 3. Furthermore the 
secondary eclipses appear to be non-existent. This can be accounted for by assuming a very 
low surface brightness for the secondary which will produce a very shallow secondary eclipse 
that could be dwarfed by the Cepheid variability. The absence of significant secondary 
eclipses could be a consequence of poor observational coverage. The near 2-day pulsational 
period (2.03 days) combined with MACHO's single-point per night coverage results in re- 
peated sampling of the same two pulsation phases during an individual eclipse. These factors, 
poor coverage and the ill-defined secondary eclipses, produce large uncertainties in the fit 
parameters and component properties in Table 11 and Figure 4. 
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In these eclipsing systems, the radial displacement of the Cepheid can be detected 
directly solely from modelling of the photometric lightcurve. This is in contrast to most 
measurements of radial amplitude which are inferred from radial velocity measurements. As 
a test, the fits of all three systems were repeated with the amplitude of the radial displacement 
held fixed at and the resulting values compared to those with radial variation included. 
For the system 6.6454.5, xl increased to 4.0, a 145% increase over the value of 1.6 listed in 
Table 9. For 81.8997.87 xl increased to 1.5, an increase of only 26% over the best fit value of 
1.1. For 78.6338.24 increased to 8.8 from 7.5, only a 17% change. For each of the systems 
the change in is found to be statistically significant for the number of degrees of freedom 
present. The relatively smaller impact of the radial ampfitude on the quafity of the fit for 
78.6338.24 could be explained by at least two factors: 

1. In this system the eclipse duration and pulsation timescales are sufficiently similar that 
they produce a degeneracy in a parameter set which includes AR. 

2. The x^ value for this system is already elevated due to the model inadequacies previ- 
ously mentioned. Their impact on the fit could easily dwarf the effects of the inclusion 
of radial amplitude. 

For 81.8997.87 the small fractional amplitude of the radial change (~6%) could make 
its effect on the lightcurve difficult to discern. 

6. SUMMARY AND FUTURE WORK 

We have analysed the available multicolor photometry for three eclipsing Cepheids. The 
principal results of this work are: 

1. MACHO Project V and R photometry for these important systems is reported. 

2. The characterization of the three systems using a self-consistent eclipsing binary lightcurve 
model which includes the effects of radius and surface brightness change due to the 
pulsating star. 

3. The interpretation of the evolutionary state of the three systems, with two of the three 
systems conforming poorly to the expectations from standard, single star evolutionary 
theory. 

4. Direct evidence of radial size change in a Cepheid variable. To date, this has been 
implied (apparently correctly!) by integrating the radial velocity curves of Cepheids. 
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Michelson interferometry is capable of measuring the angular diameters of Cepheids 
but to date the change in angular diameter, which is direct evidence for radius change, 
has only been reported for a single Cepheid (Lane et al. 2000). 

5. The prediction of future primary and secondary eclipses and their fine structure to 
facilitate follow-up observations (see below). 

6. The identification of the factors which limit the usefulness of the lightcurve model 
which has been developed. 

Follow-up photometric observations are clearly warranted for all of these systems, but 
especially for 81.8997.87 which is clearly an intermediate-mass, first-overtone Cepheid. In- 
deed, several of us have obtained multicolor photometry of the April 2001 primary minimum 
of this system using the 74-inch telescope at Mount Stromlo Observatory, Canberra, Aus- 
tralia. More densely spaced observations of this system during a future secondary minimum, 
especially longer wavelength photometry, are likely to provide significantly improved values 
for surface brightness ratios. To facilitate future observations predicted dates of primary and 
secondary minima are listed for all three systems in Tables 12-14. 

The full impact of the discovery of eclipsing Cepheids in the LMC will only be realized 
when radial velocity curves for both components in each system have been obtained. Re- 
cently, Cycle 10 HST (Hubble Space Tclccopc) observations were taken by a collaboration 
involving several coauthors on this paper and led by Edward Guinan of Villanova Univer- 
sity. For 81.8997.87, it appears that spectral observations in the far red will be necessary to 
improve the contrast of the secondary stars spectral hues relative to the Cepheid. The long 
orbital periods of these systems imply numerous HST visits over the course of approximately 
a year to map out the radial velocity curves completely. Fortunately, the relatively small 
amplitude of the radial velocity curve of the Cepheid will result in a clean separation of the 
orbital and pulsation components with a small number of visits. 

We are very grateful for the skilled support given our project by the technical staff at 
the Mount Stromlo Observatory, and in particular we would like to thank Simon Chan, Glen 
Thorpe, Susannah Sabine, and Michael McDonald for their invaluable assistance in obtaining 
the data. This work was performed under the auspices of the U.S. Department of Energy 
by the University of California, Lawrence Livermore National Laboratory under contract 
No. W-7405-Eng-48. Work performed by the Center for Particle Astrophysics personnel 
is supported in part by the Office of Science and Technology Centers of the NSF under 
cooperative agreement A-8809616. Work performed at MSSSO is supported by the Bilateral 
Science and Technology Program of the Austrahan Department of Industry, Technology and 
Regional Development. DM is supported by FONDAP. D. L. W. and D. Lepischak were 
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Table 1. V Photometry for 6.6454.5 



HJD 




V 




0-y 


Source 






(mag) 


(mag) 




2449074 , 


,95500 


14. 


.6530 





.0150 


MACHO 


2449075 , 


, 06470 


14. 


.6520 





.0150 


MACHO 


2449075 , 


,96150 


14. 


.6840 





.0150 


MACHO 


2449076. 


,99790 


14, 


.7080 





.0150 


MACHO 


2449077 . 


,11200 


14, 


.7190 





.0150 


MACHO 


2449081. 


,02760 


14, 


.7020 





.0150 


MACHO 


2449082. 


,99210 


14, 


.6240 





.0160 


MACHO 


2449083 . 


, 94080 


14, 


,6240 





.0160 


MACHO 


2449084 . 


, 97280 


14, 


.6290 





.0150 


MACHO 


2449085 . 


, 94290 


14, 


.6740 





.0160 


MACHO 


2449088 . 


, 06740 


14, 


.6020 





.0150 


MACHO 


2449088 . 


, 96770 


14, 


.6030 





.0150 


MACHO 


2449089 . 


,97160 


14, 


.6390 





.0150 


MACHO 


2449093. 


,90660 


14, 


.5990 





.0150 


MACHO 


2449098 . 


, 04000 


14, 


.6050 





.0150 


MACHO 
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Table 2. R Photometry for 6.6454.5 



HJD 




R 






Source 






(mag) 


(mag) 




2448824 . 


,17260 


14, 


.5560 





.0160 


MACHO 


2448827 . 


,18660 


14, 


.5590 





.0150 


MACHO 


2448829 . 


,14130 


14, 


.5420 





.0160 


MACHO 


2448830. 


15910 


14, 


.5260 





.0150 


MACHO 


2448833 . 


15680 


14, 


.6010 





.0150 


MACHO 


2448834 . 


,27580 


14, 


.5370 





.0150 


MACHO 


2448835 . 


, 14270 


14, 


.4810 





.0150 


MACHO 


2448836 . 


,21040 


14, 


.5340 





.0150 


MACHO 


2448837 . 


,13050 


14, 


.5520 





.0150 


MACHO 


2448843 . 


,16980 


14, 


.6000 





.0160 


MACHO 


2448844 . 


, 13490 


14, 


.5380 





.0150 


MACHO 


2448851. 


, 15840 


14, 


.5840 





.0150 


MACHO 


2448854. 


,28290 


14, 


.5080 





.0160 


MACHO 


2448855 . 


,25190 


14, 


.5290 





.0150 


MACHO 


2448857 . 


,30420 


14, 


.6080 





.0150 


MACHO 



Table 3. I Photometry for 6.6454.5 



HJD I (7/ Source 

(mag) (mag) 



2450832 , 


.78089 


14, 


.4520 


0, 


.0080 


OGLE 


2450834 , 


.74339 


14, 


.3600 


0, 


.0150 


OGLE 


2450836 , 


.68017 


14, 


.4310 


0, 


.0140 


OGLE 


2450838, 


.75705 


14, 


.3800 


0, 


.0100 


OGLE 


2450839, 


.71152 


14, 


.3690 


0, 


.0100 


OGLE 


2450840, 


.69671 


14, 


.3890 


0, 


.0070 


OGLE 


2450841, 


.77146 


14, 


.4440 


0, 


.0060 


OGLE 


2450842 , 


.69221 


14, 


.4470 


0, 


.0110 


OGLE 


2450843 , 


. 69775 


14, 


.3800 


0, 


.0100 


OGLE 


2450844 , 


.69179 


14, 


.3490 


0, 


.0070 


OGLE 


2450845 , 


.55570 


14, 


.3901 


0, 


.0079 


CTIO 


2450845 , 


.80135 


14, 


.3850 


0, 


.0080 


OGLE 


2450846 , 


.75718 


14, 


.4380 


0, 


.0070 


OGLE 


2450850, 


.73583 


14, 


.3890 


0, 


.0110 


OGLE 


2450851, 


.65660 


14, 


.4245 


0, 


.0079 


CTIO 
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Table 4. V Photometry for 78.6338.24 



HJD 




V 




0-y 


Source 






(mag) 


(mag) 




2448886 , 


, 25480 


15, 


.5390 





.0150 


MACHO 


2448888 , 


,15230 


15, 


.6790 





.0170 


MACHO 


2448889 , 


,20710 


15, 


.6940 





.0160 


MACHO 


2448893. 


, 23740 


15 


.8760 





.0170 


MACHO 


2448895 . 


.21320 


16 


.0230 





.0190 


MACHO 


2448896 . 


.18200 


16 


.1130 





.0160 


MACHO 


2448897 . 


.24090 


16 


.0870 





.0170 


MACHO 


2448902 . 


,27770 


15 


.5770 





.0150 


MACHO 


2448917, 


,16550 


15, 


.8170 





.0160 


MACHO 


2448919, 


,19520 


15, 


.4720 





.0150 


MACHO 


2448924 , 


, 20640 


15, 


.6990 





.0150 


MACHO 


2448933 , 


,25020 


16, 


.0810 





.0240 


MACHO 


2448940 , 


.15000 


15 


.6570 





.0160 


MACHO 


2448942 , 


. 17400 


15 


.7230 





.0160 


MACHO 


2448949 , 


, 13640 


16 


.1180 





.0160 


MACHO 
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Table 5. R Photometry for 78.6338.24 



HJD 




R 






Source 






(mag) 


(mag) 




2448886 , 


, 25480 


15, 


.2120 





.0150 


MACHO 


2448888 , 


,15230 


15, 


.3130 





.0170 


MACHO 


2448889 , 


,20710 


15, 


.3140 





.0160 


MACHO 


2448893. 


, 23740 


15 


.4970 





.0170 


MACHO 


2448895 . 


.21320 


15 


.6400 





.0180 


MACHO 


2448896 . 


.18200 


15 


.7300 





.0150 


MACHO 


2448897 . 


.24090 


15 


.7460 





.0170 


MACHO 


2448902 . 


,27770 


15 


.2600 





.0150 


MACHO 


2448917, 


,16550 


15, 


.5060 





.0160 


MACHO 


2448919, 


,19520 


15, 


.1810 





.0150 


MACHO 


2448924 , 


, 20640 


15, 


.3400 





.0150 


MACHO 


2448933 , 


,25020 


15, 


.7500 





.0270 


MACHO 


2448940 , 


.15000 


15 


.2980 





.0160 


MACHO 


2448942 , 


. 17400 


15 


.3340 





.0160 


MACHO 


2448949 , 


, 13640 


15 


.7660 





.0150 


MACHO 
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Table 6. V Photometry for 81.8997.87 



HJD 




V 




0-y 


Source 






(mag) 


(mag) 




2448919, 


, 24260 


17, 


.1530 





.0220 


MACHO 


2448920 , 


, 04760 


17, 


.0300 





.0230 


MACHO 


2448921 , 


, 04080 


17, 


.1880 





.0220 


MACHO 


2448924 . 


,25670 


16, 


.9930 





.0230 


MACHO 


2448927 . 


.99640 


16, 


.9700 





.0270 


MACHO 


2448929 . 


.03000 


17, 


.1960 





.0200 


MACHO 


2448929 . 


.97740 


16, 


.9830 





.0240 


MACHO 


2448930 . 


.97620 


17, 


.1670 





.0360 


MACHO 


2448931 , 


,99100 


17, 


.0060 





.0240 


MACHO 


2448932 , 


, 16440 


16, 


.9770 





.0240 


MACHO 


2448934 , 


, 13420 


16, 


.9560 





.0220 


MACHO 


2448936 , 


, 13460 


16, 


.9940 





.0420 


MACHO 


2448938 . 


, 18980 


16, 


.9860 





.0300 


MACHO 


2448939 . 


.01450 


17, 


. 1340 





.0280 


MACHO 


2448940 , 


.03150 


17, 


.0110 





.0240 


MACHO 
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Table 7. R Photometry for 81.8997.87 



HJD 


R 


(Tr Source 




(mag) 


(mag) 



2448919, 


, 24260 


16, 


.4420 


0, 


.0170 


MACHO 


2448920 , 


, 04760 


16, 


.3470 


0, 


.0180 


MACHO 


2448921 , 


, 04080 


16, 


.4710 


0, 


.0170 


MACHO 


2448924 . 


,25670 


16, 


.3450 


0, 


.0180 


MACHO 


2448927 . 


.99640 


16, 


.3040 


0, 


.0190 


MACHO 


2448929 . 


.03000 


16, 


.4600 


0, 


.0170 


MACHO 


2448929 . 


.97740 


16, 


.3200 


0, 


.0180 


MACHO 


2448930 . 


.97620 


16, 


.4390 


0, 


.0210 


MACHO 


2448931 , 


,99100 


16, 


.3390 


0, 


.0180 


MACHO 


2448932 , 


, 16440 


16, 


.3170 


0, 


.0180 


MACHO 


2448934 , 


, 13420 


16, 


.2960 


0, 


.0170 


MACHO 


2448936 , 


, 13460 


16, 


.3140 


0, 


.0240 


MACHO 


2448938 . 


, 18980 


16, 


.3100 


0, 


.0210 


MACHO 


2448939 . 


.01450 


16, 


.4250 


0, 


.0190 


MACHO 


2448940 , 


.03150 


16, 


.3440 


0, 


.0180 


MACHO 



Table 8. I Photometry for 81.8997.87 



HJD I (7/ Source 

(mag) (mag) 



2450739 , 


.86067 


15, 


.6720 


0, 


.0100 


OGLE 


2450744 , 


.76901 


15, 


.7840 


0, 


.0130 


OGLE 


2450745 , 


.85362 


15, 


.6680 


0, 


.0130 


OGLE 


2450746, 


.83950 


15, 


.7950 


0, 


.0180 


OGLE 


2450747 , 


.73495 


15, 


.6580 


0, 


.0130 


OGLE 


2450750, 


.79988 


15, 


.8020 


0, 


.0130 


OGLE 


2450751, 


.80069 


15, 


.6610 


0, 


.0200 


OGLE 


2450752, 


.85852 


15, 


,7830 


0, 


.0110 


OGLE 


2450755 , 


.78274 


15, 


.6530 


0, 


.0100 


OGLE 


2450759 , 


.70861 


15, 


.6680 


0, 


.0150 


OGLE 


2450761 , 


.78982 


15, 


.6860 


0, 


.0130 


OGLE 


2450766 , 


.83954 


15, 


.7490 


0, 


.0150 


OGLE 


2450773 , 


.79513 


15, 


.7050 


0, 


.0140 


OGLE 


2450776 , 


.82033 


15, 


.7410 


0, 


.0140 


OGLE 


2450778 , 


.85342 


15, 


.7100 


0, 


.0100 


OGLE 
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Table 9. Best fit parameters for 6.6454.5. Meanings of individual parameters and units 

are explained in the text. 







1.6 






Porb 


id) 


397.142 ± 0.005 






in 




86.71 ± 0.04 






Companion 


(primary) 


Variable 


(secondary) 


R 




0.0522 ± 0.0004 


Rfnin 


0.051 ± 0.001 


Jv 
Jr 




1.061 ± 0.002 


<*> 


0.748 ± 0.005 


Xv 




0.46 ± 0.17 


Xv 


0.52 ± 0.85 


Xr 




0.48 ± 0.20 


Xr 


0.42 ± 0.72 


Xl 




0.45 ± 0.18 


Xl 


0.56 ± 0.55 








Pceph {d^ 


4.97371 ± 0.00002 








^Ramp 


0.0166 ± 0.0002 








^Rshift (d) 


-0.17 ± 0.01 


Intensity-weighted Mean Magnitudes and Colours 


V 




14.82 ± 0.04 


<V> 


16.7±0.1 


R 




14.80 ± 0.04 


<R> 


16.3± 0.1 


I 




14.72 ± 0.04 


<I> 


15.89± 0.09 


V - 


R 


0.018 ± 0.002 


<V-R> 


0.452 ± 0.007 


V - 


I 


0.102± 0.003 


<V-I > 


0.81± 0.01 








<Wr> 


14.4±0.1 
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Table 10. Best fit parameters for 78.6338.24. Parameter definitions and units are the 
same as in Table 9 with the exception of Ai and A2 which are defined by equation 11. 





7.5 




Porb {d) 


419.718 ± 0.008 






86.94 ± 0.02 




Companion 


(primary) Variable 


(secondary) 



R 

Jv 
Jr 

Xv 
Xr 



0.0403 ± 0.0002 
1.146 ± 0.001 
0.46 ± 0.17 
0.48 ± 0.20 



Rfnin 

Xv 
Xr 

Pceph (^) 
^1 

A2 id) 

^Ramp 

ARshift (d) 



0.0717 ± 0.0004 
0.9118 ± 0.001 
0.52 ± 0.85 
0.42 ± 0.72 
17.68586 ± 0.0003 
2.265 ± 0.007 
-448 ± 2 
0.03106 ± 0.0002 
-6.07 ± 0.01 



Intensity-weighted Mean Magnitudes and Colours 



V 


16.55 ± 0.02 


<V> 


16.20±0.03 


R 


16.35 ± 0.02 


< R> 


15.75± 0.03 


V-R 


0.20 ± 0.03 


<V-R> 


0.447 ± 0.001 






<Wr> 


13.98±0.03 
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Table 11. Best-fit parameters for 81.8997.87. Parameters are tfie same as in Table 9. Tlie 
large uncertainties in the limb-darkening coefficients of the secondary indicate that while 
these parameters are allowed to vary, they are essentially unmodified and unconstrained by 
the lightcurve. The tabulated xx are the original values and the errors are included for 

consistency. 





1.2 






Porb {d) 


800.5 ± 0.1 






in 


87.0 ± 0.2 






Variable 


(primary) 


Companion 


(secondary) 


Rmin 


0.029± 0.001 


R 


0.047 ± 0.005 




1.078 ± 0.009 


Jv 
Jr 


0.54 ± 0.05 


Xy 


0.50 ± 0.94 


Xv 


0.50 ± 6.06 


Xr 


0.50 ± 0.91 


Xr 


0.50 ± 4.52 


Xi 


0.50 ± 0.98 


Xi 


0.50 ± 3.02 


Pceph (c^) 


2.035321 ± 0.000009 






^ Ramp 


0.0017 ± 0.0002 






^Rshift (d) 


0.21 ± 0.01 






Intensity-weighted Mean Magnitudes and Colours 


< > 


17.2 ± 0.2 


V 


20±1 


< R> 


16.6 ± 0.2 


R 


19.2± 0.9 


<I> 


16.0 ± 0.2 


I 


17.9± 0.8 


<V-R> 


0.602 ± 0.009 


V-R 


1.3± 0.1 


<V-I > 


1.137± 0.007 


V -I 


2.6± 0.2 


<Wr> 


14.1±0.2 
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Table 12. Predicted dates of future eclipses for 6.6454.5. 



Primary Eclipse Secondary Eclipse 

JD UT JD UT 



2452454.78 
2452851.95 
2453249.12 
2453646.28 

2454043.45 
2454440.62 
2454837.78 
2455234.95 
2455632.12 



2002 Jun 28 18.79 

2003 Jul 30 22.80 

2004 Aug 31 2.80 

2005 Oct 2 6.81 

2006 Nov 3 10.81 

2007 Dec 5 14.82 

2009 Jan 5 18.82 

2010 Feb 6 22.83 

2011 Mar 11 2.83 



2452653.37 
2453050.53 
2453447.70 
2453844.87 

2454242.03 
2454639.20 
2455036.37 
2455433.53 
2455830.70 



2003 Jan 13 8.80 

2004 Feb 14 12.80 

2005 Mar 17 16.81 

2006 Apr 18 20.81 

2007 May 21 0.82 

2008 Jun 21 4.82 

2009 Jul 23 8.83 

2010 Aug 24 12.83 

2011 Sep 25 16.84 



Table 13. Predicted dates of future eclipses for 78.6338.24. 



Primary Eclipse Secondary Eclipse 

JD UT JD UT 



2452379.74 
2452799.20 
2453218.67 
2453638.14 
2454057.60 
2454477.07 
2454896.53 
2455316.00 
2455735.46 



2002 Apr 14 17.70 

2003 Jun 8 4.88 

2004 Jul 31 16.07 

2005 Sep 24 3.25 

2006 Nov 17 14.43 

2008 Jan 11 1.61 

2009 Mar 5 12.79 

2010 Apr 28 23.98 

2011 Jun 22 11.16 



2452589.47 
2453008.94 
2453428.40 
2453847.87 
2454267.33 
2454686.80 
2455106.27 
2455525.73 
2455945.20 



2002 Nov 10 11.29 

2004 Jan 3 22.48 

2005 Feb 26 9.66 

2006 Apr 21 20.84 

2007 Jun 15 8.02 

2008 Aug 7 19.20 

2009 Oct 1 6.39 

2010 Nov 24 17.57 
2012 Jan 18 4.75 
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Table 14. Predicted dates of future eclipses for 81.8997.87. 



Primary Eclipse Secondary Eclipse 

JD UT JD UT 



2452010.81 
2452811.73 
2453612.64 
2454413.56 
2455214.47 
2456015.39 
2456816.30 
2457617.22 
2458418.13 
2459219.05 



2001 Apr 10 19.47 
2003 Jun 20 17.44 
2005 Aug 29 15.40 
2007 Nov 8 13.37 
2010 Jan 17 11.33 
2012 Mar 28 9.30 
2014 Jun 7 7.26 
2016 Aug 16 5.23 
2018 Oct 26 3.20 
2021 Jan 4 1.16 



2452411.27 
2453212.18 
2454013.10 
2454814.01 
2455614.93 
2456415.85 
2457216.76 
2458017.68 
2458818.59 
2459619.51 



2002 May 16 6.45 
2004 Jul 25 4.42 
2006 Oct 4 2.38 
2008 Dec 13 0.35 
2011 Feb 21 22.32 
2013 May 2 20.28 
2015 Jul 12 18.25 
2017 Sep 20 16.21 
2019 Nov 30 14.18 
2022 Feb 8 12.14 
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Fig. 1. — Primary eclipses of 6.6454.4 in V with curve of best fit. Upper panels show 
residuals in magnitudes. Filled boxes indicate observations from the MACHO project, open 
boxes indicate observations from the OGLE project, x's are follow-up observations taken at 
CTIO. 
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Fig. 2. — Primary eclipses of 78.6338.24 in V with curve of best fit. Upper panels show 
residuals in magnitudes. All points are from the MACHO project database. 
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Fig. 3. — Primary eclipses of 81.8997.87 in V with curve of best fit. Upper panels show 
residuals in magnitudes. Filled boxes indicate observations from the MACHO project, open 
boxes observations from the OGLE project. 
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Fig. 4. — Wr vs. logioP diagram for 1766 MACHO Cepheids. The brighter sequence at a 
given period are Cepheids pulsating in the first overtone and the sequence extending to longer 
periods are fundamental mode pulsators. Stars in the lower right arc Type II Cepheids. The 
locations of the three Cepheids studied here are indicated. Their magnitudes and colors are 
from the best-fit parameters and thus have had the companion flux removed. 
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Fig. 5. — Systems 6.6454.5 and 81.8997.87 compared with theoretical isochrones for Y = 
0.25, Z = 0.008 stars (BerteUi et al. 1994) ranging from log(age)=9.1 to log(age)=7.9 [years]. 
Also shown are the theoretical fundamental (solid) and overtone (dashed) instability strips 
of Chiosi, Wood, & Capitanio (1993). Cepheids are shown as solid squares, their companions 
as solid triangles. Dashed lines connect system members. 



-35- 




Fig. 6. — Systems 6.6454.5 and 78.6338.24 compared with theoretical isochrones for Y = 
0.230, Z = 0.0004 stars (Fagotto, Bressan, BerteUi, & Chiosi 1994). Also shown are the 
theoretical fundamental (sohd) and overtone (dashed) instabihty strips of Chiosi, Wood, & 
Capitanio (1993). 



